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Abstract 

We investigate the stationary state of a model system evolving according to a modified focusing 
truncated nonlinear Schrodinger equation (NLSE) used to describe the envelope of Langmuir waves 
in a plasma. We restrict the system to have a finite number of normal modes each of which is in 
contact with a Langevin heat bath at temperature T. Arbitrarily large realizations of the field are 
prevented by restricting each mode to a maximum amplitude. We consider a simple modeling of 
wave-breaking in which each mode is set equal to zero when it reaches its maximum amplitude. 
Without wave-breaking the stationary state is given by a Gibbs measure. With wave-breaking the 
system attains a nonequilibrium stationary state which is the unique invariant measure of the time 
evolution. A mean field analysis shows that the system exhibits a transition from a regime of low 
field values at small |A|, to a regime of higher field values at large |A|, where A < specifies the 
strength of the nonlinearity in the focusing case. Field values at large |A| are significantly smaller 
with wave-breaking than without wave-breaking. 
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I. INTRODUCTION 



The statistical mechanics of continuum fields is a problem of both intrinsic and applied 
interest. It is often at the level of the coarse-grained continuum description provided by the 
collective modes of the system, rather than at the underlying level of the atomic degrees 
of freedom, where certain types of physical phenomena are most clearly manifest. Thus, 
turbulence in fluids is best described by statistical properties of hydrodynamic fields evolving 
according to the Navier-Stokes equations.. 

Another example, relevant for plasmas, is the case of a complex field the dynamics 
of which is governed by the nonlinear Schrodinger equation: 

id t <f) + A0 = A|0| p " 2 0, 

(1) 

t > 0, x eAcR D , and <p{x,t = 0) = (f) (x), 

where A G K is a parameter specifying the strength of the nonlinearity, p > 2 is an integer, 
and A is a compact domain in M. D which we shall generally take to be a torus. Equation 
(1) with A < (the focusing case) and p = 4 has been used to describe some of the physics 
associated with the slowly varying envelope of Langmuir waves in a plasma, such as wave 
collapse in the subsonic regime [1], as well as other physical phenomena [2]. For A < the 
Hamiltonian generating the dynamics, 

^(0) = ^Q|V0| 2 + ^|0| p ) d D x, (2) 

is not bounded below. Consequently there does not exist a Gibbs measure for the fields 
evolving according to (1). One way of overcoming this problem is by putting bounds on the 
conserved "mass", 

M<P) = \J A \<P\ 2 d D x. (3) 

This is done in [3] for D — 1, see also [4], where a well-defined measure is constructed for 
parameter ranges where there is no wave collapse [1]. This approach is difficult to extend 
to D > 1 where one is faced with ultraviolet problems, see [5] for the defocusing case, and 
is not suitable to cases where wave collapse can occur. 

In this paper we investigate the properties of a field evolving according to a truncated 
version of (1) with added stochastic terms as well as an upper bound on the amplitude of 
the Fourier components of 0. The stochastic Langevin forces and the focusing nonlinearity 
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want to drive the field to arbitrarily large values, which is prevented by the bound on the 
amplitude. We also consider a mechanism which resets the amplitude of a Fourier mode 
to zero when it reaches its maximum. This mechanism imitates, in a simplified way, the 
wave-breaking of Langmuir waves in a plasma [6] [7], which we now describe. 

The wave-breaking of a Langmuir wave is the violent dissipation of its energy when 
its amplitude is high enough for the bulk of the thermal electrons to get trapped in the 
electrostatic potential of the wave [6] [7]. The larger the wave number k the smaller the 
maximum amplitude the wave can reach. For k > A^ 1 , where \d is the Debye screening 
length, wave-breaking turns into strong Landau damping so the amplitude does not grow. 
Wave-breaking introduces both a natural ultraviolet cut-off, k max — A^, 1 , and a natural 
upper bound for the amplitude of each Fourier mode of 0. Our modeling of wave-breaking is 
a discontinuous dissipative process whereby each Fourier component of (f> is set equal to zero 
when it reaches a given amplitude. In the absence of wave-breaking the stationary measure 
is Gibbsean at the temperature imposed by the Langevin heat bath. With wave-breaking 
the stationary state is not an equilibrium one. We show however that it exists and is unique. 

The outline of the paper is as follows. In Section II we specify the stochastic dynamics 
of our model system. This dynamics acting on a finite number of bounded degrees of freedom 
(the Fourier modes) ensure the existence of a unique stationary state which is approached 
exponentially fast. This is proven in Section III. A mean field theory with and without 
wave-breaking is developed in Section IV. Consequences of an infrared divergence of the 
maximum amplitude in the mean field theory with wave-breaking are briefly addressed in 
Section IV D. Finally, the relevance of the mean field theory and some possible improvements 
of our wave-breaking prescription are discussed in Section V. 

II. MODEL AND DEFINITIONS 

Let A be a D-dimensional torus of length L and volume V = L D . From now until Section 
IV C, we will take L — 1 without loss of generality. Let r\ c > be an ultraviolet cut-off [in 
plasma physics, r] c ~ (2n\ D )~ 1 ], i.e. consider a field 4> of the form 




(4) 



n\\<rj c 
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with k n = 2nn, n G Z D , and ||n|| = (J2?=i n l Y^ 2 - The truncated Hamiltonian and mass are 
given respectively by 

v 

d D x, (5) 
Af(a) = N{<t>) = \ Kl 2 - (6) 

||n||<7? c 

Note that T-L(a) and Af(a) are different from their untruncated counterparts "H(a) and J\f(a) 
defined in (2) and (3). In the absence of a bound on |a n (t)|, the time evolution of the a„ 
is assumed to be given by the truncated nonlinear Schrodinger equation (1) combined with 
stochastic Langevin forces: 

{a nR = dU(a)/da nI - v n d\pAf(a) + U{a)]/da nR + v /2u n (3~ 1 L nR , 
a„i = -dU(a)/da nR - v n d[(jj\f(a) + V.(a)]/da nI + v / 2v n (3- l L nI , 

where a nR = Re(a n ), a n i = Im(a n ), v n is the linear damping of a n , (3 = (kbT)~ 1 , T is the 
temperature, \x is the chemical potential, and the L n are independent Gaussian white noises 
with (L nR (t)L nI (t')) = and (L nR (t) L nR (t')) = (L nI (t)L nI (t')} = 5(t-t'). 

The motivations for (7) are twofold : first, in the limit of no contact with the thermostat 
(i.e. v n = for all n), it reduces to a truncated version of (1) with of the form (4) ; secondly, 
in the well defined defocusing case, A > 0, with an ultraviolet cut-off on n, the stationary 
probability distribution of the a n is the Gibbs distribution P eq (a) oc e~^ H ^ + ^ r ^\ where 
(3 and \i are specified by the Langevin terms [8]. These can be thought of as representing 
the interaction between the waves and the particle degrees of freedom, not included in (5). 
Even with this truncation the Hamiltonian "H(a) in (5) is still not bounded below for the 
case A < 0, the focusing case. Hence the dynamics in (7) will not lead to a stationary state: 
the reservoirs would just keep extracting energy from the waves, causing "H(a) to go towards 
— oo. To prevent this we let a n be a monotone decreasing function of such that a n = 
for all n with ||n|| > i] c and a n > otherwise. Now, we assume that, as soon as \a n \ = a n , 
it is either reflected, i.e. a n is set equal to — a n , or instantaneously set equal to a n = 
("wave-breaking" prescription) from where it resumes its time evolution according to (7). 
These prescriptions are equivalent to imposing either a reflecting or an absorbing boundary 
condition at \a n \ = a n on the Fokker-Planck equation associated with (7). In the latter case 
we compensate exactly for the absorption with sources at a n = 0. 
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H{a) =%{<!>) = \ J2 k 2 n K\ 2 + - [ 
A ,,"777 V J A 



^2 a n exp(ik n 



■ x) 



\n\\<Vc 



\n\\<Vc 



The probability distribution P(a, i) of the above system satisfies the modified Fokker- 
Planck equation 

d t P(a,t) - LP(a,t) = e ^ a n ({a m ^ n },t) S(a nR )5(a nI ), (8) 

IMI<»?c 

with e = in the reflecting case and e = 1 in the wave-breaking case. The Fokker-Planck 
operator L is given by 

LP(a,t) = - V„-J„(a,t), (9) 

||n||<j? c 

with the probability current 

J n (a,t) = {V n xH(a)-is n V n [H(a)+tiN(a)]}P(a,t) 

- z/jr^PM). (io) 

The sources a n ({a m ^ n },t) are given by 

o n ({a m ^n},t) = -v n ^ x \ e n -V n P(a,t)ds. (11) 

J\a n \=a n 

Here we have written V„ = (d/da n R, d/da n i), V„ x H(a) = [dH{a)/da n i,—dH(a)/da n R], 
e n = (a n R/\a n \, a n i/\a n \) and ds = \Jdo? nR + da^j. The boundary conditions on the so- 
lution to (8) are e n ■ J n (a, 0l|a n |=a n , \a m ^ n \<a m = for all t > in the reflecting case, 
and P(cL,t)\\a n \=a n ,\a m ^ n \<a m = for all t > in the wave-breaking case. We also have 
P(a, t) | \a„\>a„ = f° r an n an d all t > in both cases. 



III. EXPONENTIAL APPROACH TO A UNIQUE STATIONARY SOLUTION 

In this section we show that the system with time evolution described by (8) has a 
unique stationary state, P s t(a), and that every initial probability distribution tends to 
P s t{a) exponentially fast. For the reflecting case P s t{a) is the equilibrium Gibbs state, 
P s t(a) oc e-Pi' H ( a )+t MN '( a )} for \a n \ < a n , while for the wave-breaking prescription P s t(a) will 
be a nonequilibrium stationary state which is unknown in general. Here we give the proof 
for the wave-breaking prescription only as the equilibrium case is well known [8] . 

The Langevin equations (7) with the wave-breaking prescription defines a Markov 
process on the compact domain T = supp ITI| |<»j c l{KI<a n }- Let P(a,t;b), with a,b G T, 
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denote the solution to (8) with P(a, 0; b) = 5 (a — b) (with obvious notation), and write 
P(A, t;b) = J aeA P(a,t;b) d N a. Doeblin's condition [9] says that if 

p(t) = sup \P(A,t;b)-P(A,t;b')\<l, (12) 
6,6'er, Acr 

for some t > 0, then there is a unique invariant measure with exponentially fast convergence 
to it, starting from any point in T. This follows from the fact that, using the Chapman- 
Kolmogorov equation for Markov processes, we have p(t + s) < p(t)p(s) [9]. This implies in 
turn that if p(t) satisfies (12) for some t > then p(t) — > exponentially as t — >■ +oo. In 
order for (12) to be violated it is necessary that ther exists a set A and a configuration b 
for which P(A,t; b) = for all t > and of course P(A c ,t; b) = 1. Hence to prove (12) it is 
clearly sufficient to show that 

inf P(a,t;b) > c> (13) 
ber 

for every b G T and some t — to > 0, and all a in some compact set ifcTof strictly positive 
Lebesgue measure. When this is true then clearly c| A n if | < -P(^4, £; fe) < 1 — c| A c D if | and 
p(t) < 1 -c|if| < 1. 

To prove (13) for our system we simply note that for a and b away from dT, the 
boundary of T, i.e. |a n | < a n , we always have P(a, 1;6) > Pdiff(a,l;b) where P di ff(a,t;b) 
is the transition probability density for the evolution described by (8) without sources (i.e. 
a n — for all n, and absorbing boundaries at \a n \ = a n ). This evolution would not conserve 
mass but Pdiff(a, 1; b) > c > for a, b away from dT. Letting now b be arbitrary close to 
dT one has P(a, 1; 6) — >■ P(a, 1; 6 pr ) where b pr is the projection of b according to the wave- 
breaking prescription, which is away from dT, and the proof is complete. (We are indebted 
to S. Varadhan for this streamlined argument). 

As an illustration, Figure 1 shows the stationary distribution of /3 1 / 2 |a n | obtained from 
a numerical solution of (7) for a three-mode system (n = 0, ±1) with D — 1, p — 4, 
= —5, v n = 10, (3 l / 2 atQ = 1, and (3 l l 2 a±i = 1/2. For each mode the stationary 
distribution is computed from the relative frequency of /3 1//2 |a n | G [5m 1CT 3 , 5(m + 1) 1CT 3 [, 
with m an integer, over a long enough period of time. The time at which we start the 
sampling is chosen large enough so that the results do not depend on its value. Figure 2 is 
the counterpart of Figure 1 for reflecting boundary conditions. Analytical results obtained 
from the corresponding Gibbs distribution are also shown (plain lines). In this particular 
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(b) 




0.0 0.2 0.4 0.6 0.8 1.0 



/3 l/2 K\ 

FIG. 1: Stationary distribution of (a) /3 1//2 |ao| and (b) /3 1 / 2 1 a-i | obtained from a numerical solution 
of (7) for a three-mode system (n = 0, ±1) with D = 1, p = 4, A/3 -1 = -5, v n = 10, Z? 1 / 2 ^ = 1, 
and /3 1//2 ay-i = 1/2. Tiie distribution for /3 1//2 |a_i| is similar to (b). 

case the effects of the wave-breaking prescription are most clearly seen in the distribution 
of | do | - For reflecting boundary conditions (Fig. 2), the distribution of |a | is maximum 
at the boundary, which means that this mode tends to stay in the vicinity of its maximum 
amplitude for a relatively long period of time. On the other hand (Fig. 1), the wave-breaking 
prescription prevents the mode from staying near its maximum amplitude very long. Instead, 
it is quickly reset to \a \ =0, which results in a dramatic reduction of the distribution near 
the boundary and the appearance of a narrow peak at |ao| =0. 

IV. MEAN FIELD THEORY (p = 4) 

We turn now to a mean field approximation similar to that considered by various au- 
thors in various contexts (see e.g. [10]). It corresponds to replacing ^\<fi(x)\ p in (2) with 
l,W(x) p l 2 ~ l \(j)(x)\ 2 where W(x) is to be determined self-consistently in a way to be dis- 
cussed below. In the following, we will take p — 4, assume that W(x) = W is independent 
of x, and consider both the reflecting (e = 0) and wave-breaking (e = 1) cases. 
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0.0 0.2 0.4 0.6 0.8 1.0 



/3 m K\ 

FIG. 2: Stationary distribution of /3 1//2 1 a.^ | obtained from a numerical solution of (7) with reflecting 
boundary conditions, and analytical results obtained from the corresponding Gibbs distribution 
(plain lines). Parameters and notation are the same as in Fig. 1. 



A. General results 



Our mean field theory is defined by replacing / H((f>) with 

n MF {<t>) = jf Q|V0| 2 + ^|0| 2 ) d D x. 



(14) 



In this approximation, the Fourier components of <fi are decoupled and evolve according to 
(7) in which H(a) is replaced with 

H M F(a) = Hmf(4>) = \ J2 {kl + ^W)\a n \ 2 . 



\n\\<Vc 



According to (8), the mean field stationary measure of a n is the solution to 



= 
+ 



d 



d 



da nI 



dn M F(a) _ 9Hmf(o) +liN(a) . p 
oa nI oa nR 

dW,MF(a) &Hmf{o) + fjJ\f(a 



da n R 
d 2 d 2 



+ v n - 



P(a n ) 



da l R da nl 



da nI 

2 i P(a>n) + £o- n 5(a nR )5(a nI ), 



with the boundary conditions 



e„ • Jn(a n )\\a n \= an = 0, fore = 0, 
P{a n )\\a n \=a n = 0, fore = l, 



(15) 



(16) 



(17) 
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and where 



o n = -v n fl 1 e n ■ V„P(a„) ds. 

J\a„\=a n 



(18) 



Let FmfK) denote this stationary measure. One has, 



PMF\P"ni 



l -\a n \<a„ 

zjwj 

1 



x 



exp 
1 



s\a n \/a„ r 



cxp 



e/3 



K + XW + 



j a 



2 2 

r 

n 



dr, 



(19) 



where Z n (W) is a normalization constant such that f, , . -Pmf(gO da n nda n i = 1, and the 
integral over r is equal to 1 when e — 0. Before writing the expression of Z n (W) it is 
convenient to introduce the function 



M*) = f (*£ + As + |). 

After some algebra one finds that Z n (W) is given by 



7T 



[l-exp(-aXW)], 



in the reflecting case (e = 0), and by 

^ [Ei^X W) - ln(a^ n(W ,)) _ 7 ] f or > , 

2\h^w)\ i E ~i( a n\ h n(W)\) + \n(a 2 n \h n (W)\) + 7] for ^(W) < 0, 



(20) 



(21) 



Z n (W) = 



(22) 



where 7 is the Euler constant, in the wave-breaking case (5 = 1). 



B. Self-consistency and nonuniqueness 

There are various self-consistent prescriptions for finding W and it is a matter of judge- 
ment which gives the best approximation. Two "natural" ways to impose self-consistency 
are: (i) require that the average energies agree, which yields 

\w^~\\^) w = l -(W)w, 
z p 

where (-)w denotes the average w.r.t. the stationary mean field measure; and (ii) require 
that in the dynamics specified by (1), 

w^ 2 -v = (i0r 2 w- 

In the equilibrium case corresponding to e = 0, a third "natural" prescription is: (iii) require 
that W minimize the free energy (24). 
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For the case p = 4, A > 0, a n ,\\n\\<r lc = +00, and a Gaussian stationary measure, these 
prescriptions give the self-consistency equation W = q(\4>\ 2 )w with q — 1 for (i) and (ii), 
and q = 2 for (hi). 

We shall use the same self-consistency equation in our model also for A < and finite 



n\\<r) c 



w with (lanl 2 )^ — J^ <an \a n \ 2 PMF(cin) da n Rda n i. 



«n,||n||<^, giving W = qY,\ 
We shall leave q as a free parameter from now on as the results depend on the choice of q 
only via a rescaling of A. Using (19), (21), and (22), this yields the self-consistency equation 



W = q 



\n\\<Vc 



K(w) 



1 - 



exp((5 - l)a 2 n h n {W)) 



(23) 



¥Z n (W) 

We now consider the case where (23) has more than one solution. We then need to 
determine which one is to be picked. In the reflecting case this can be done by choosing the 
solution with the smallest free energy 



E ■ 

|n||<r? c 



+ 



Pmf(o>) In Pmf(o) da, 



(24) 



w 



which, up to a term independent of W, is equal to 

4 

1 



A 



E a - 

\n\\<r) c 



) -\w E <ki 2 >J-I E 

I w Nlofc / p |M|o,c 



(25) 



As mentioned above, one could alternatively define W, in the reflecting the value 

minimizing (25). It will be shown later that, in the large system limit, this definition of W 
does not give significantly different results from our definition with q = 2 (i.e. non Gaussian 
corrections are negligible in this limit). 

In the wave-breaking case, the stationary state is a nonequilibrium state and there is 
no known counterpart of (24) that would make it possible to choose among the solutions to 
(23) in a simple way. Under these circumstances, we will give all the solutions and base our 
discussion on the well-defined limits of a small or large |A| where (23) has a single solution. 



C. Determination of W for A < in the large system limit 

We now solve (23) with A < in the case of a large system for both e = and e — 1. 
We show that the mean field theory predicts a transition from a regime of low W at small 
|A| to a regime of higher W at large |A|, and that the value of W at large |A| is significantly 
smaller for e — 1 than e — 0. 
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From now on we no longer take L — 1. This will make V appear explicitely in the self- 
consistency equation (23). Normalizing the a n so that the right-hand side of (4) is multiplied 
by V -1 / 2 , it can be seen that the right-hand side of (23) must be multiplied by V^ 1 , and a n 
must be replaced with V 1 ^ 2 a n on the right-hand side of (21), (22), and (23). Note also that 
in the expression (20) for h n (x), k n is now given by k n = 27in/L, n G IP '. In the following 
we will take a n = a(\\k n \\) with a(k) the max of zero and a linearly decreasing function of 
k= \\k\\. 

We consider the large system limit in which V is large enough so that discrete sums 
over n G IP can be replaced by integrals over k G MP according to 



V ^ (2ttV 



Fix a > 0, k max > 0, and let 



, , Q!o(l k/k max ) if k < k max , 
a{k) = < (27) 



Make the change of variables 



and define 

S D 



if k > k ma n 



x = /i/2- \X\W, 
m = f/,/2, 



(28) 



(2tt) d J x + k 2 



(3Va{kf 
Z e (x,k,l3V) 



k°- l dk, (29) 



in which Sd denotes the surface of the unit D-sphere, 




Z (x, k, PV) = |exp[^\/a(A;) 2 (x + A; 2 )] - 1 J> , (30) 



^{^[iPVaikfix + k^-XnllpVaikfix + k 2 )}-^ for x > -k\ 

J ^{E 1 [^Va(k) 2 \x + k 2 \] + \n[^Va(k) 2 \x + k 2 \}+ 1 } for x < -k 2 . 

(31) 

Using (26), (28) and (29), one can rewrite (23) as 

5 | A i 

m — x — — —T e (x, f3V) (x<m), (32) 

the solution to which is of the form x = x £ (m, g|A|/3 _1 , (5V). From this result, (28), and the 
self-consistency equation W = q(\(fi\ 2 )w we finally find 

PW') = fMMnV) S "-y^ . (33) 
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FIG. 3: fi{\4>\ 2 ) as a function of q\\\P~ l [Eq. (33)] in the reflecting boundary case (e = 0) for 
D = 3, n = 1 (i.e. m = 1/2), «o = 1> k max = 1, and j3V = 10 3 . The solid line corresponds to the 
solution with the smallest free energy (25). 

Figure 3 shows a typical example of /3(\(f)\ 2 ) as a function of g|A|/3 _1 in the reflecting 
case (e — 0) for D — 3, // = 1 (i.e. m = 1/2) corresponding to the value for Langmuir 
waves, ao = 1, k max = 1, and /3V = 10 3 . When (32) has more than one solution, we pick 
the one with the smallest free energy (25) (solid line). One observes a sharp transition at 
g|A|/9 _1 = 0.69 ± 0.01 from a low (|0| 2 ) at small |A| to a higher (|0| 2 ) at large |A|. Similar 
results obtained for different values of /3V, all the other parameters being fixed, show that 
for small |A| one has the usual scaling (|0| 2 ) ~ /3' 1 . In this regime, the modes are stable 
and localized near \a(k)\ = where they are not significantly affected by the bound on 
their amplitude. In the opposite case of large |A|, the modes are unstable and localized near 
their maximum \a(k)\ ~ \fV . This is a saturated regime, independent of |A|, and one has 
the simple scaling (|0| 2 ) ~ V. The value of A at which the transition occurs scales like 

|A| ~ (qvy 1 . 

The curve for /3(|0| 2 ) as a function of 2|A|/3 _1 for W minimizing (25) is indistinguishable 
from the solid line in Fig. 3 with q = 2, which means that minimizing (25) gives practically 
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the same result as W = 2(\(f)\ 2 } w . This is a consequence of the large system limit in which 

1 





4 




> -1 


\\n\\<Vc 


' w 



J {Hk)\ 2 )wd D k 



+ 



1 



~ 2 



V(2irY 



1 



n / [(\a(k)\*) w -2(\a(k)\Yw] d D k 



(2tt)' 



„ I (\a(k)\ 2 ) w d D k 



Using this expression in (25) and minimizing w.r.t. W one obtains, 

2 



W 



(27T) 



^ J(\a(k)\ 2 ) w d D k = 2(\cf ) \'< 



)w- 




Figure 4 is the counterpart of Figure 1 for the wave-breaking case (e — 1) with the same 
parameters. The behavior of the curve at low and large |A|, where there is only one solution, 
clearly suggests the existence of a transition, similar to the one observed in the reflecting 
case, at some intermediate |A|. The most important difference is the strong reduction of 
(|0| 2 ) in the large |A| regime. Because of the wave-breaking prescription, unstable modes 
cannot remain localized near their maximum very long, as they do in the reflecting case. 
Instead they are rapidly reinjected at \a(k)\ = 0, which causes the observed reduction of 
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(|0| 2 ) compared with Fig. 1. It can also be seen that (|0| 2 ) depends on |A| in this regime, 
which is thus not saturated. This is due to the fact that an unstable mode is more likely to 
hit its boundary, and be reset to \a(k)\ = 0, at larger |A|. The average time it spends near 
its maximum is reduced, which explains the observed slow decreasing of (|0| 2 ) at large |A|. 

It must be emphasized that the phase transitions reported here are a characteristic of 
the mean field theory which does not exist in the exact theory for a finite V. The result of 
Section III rules out any possibility of coexistence of distinct stationary solutions to (8) at 
any specified A. 



D. Mean field theory with an infrared diverging a(k) and e = 1 in the large system 
limit 

For Langmuir waves in plasmas, one has both wave-breaking (e — 1) and an infrared 
diverging a(k) behaving like 1/k for k — > 0. Careful estimates give [6] [7] 

a{kf oc ^ , (34) 

for k < l/y/3 and a(k) = for k > l/y/S. The simplest way of getting meaningful results 
in this setting consists in getting rid of the unbounded k = mode by forcing a (t) = for 
all t > 0. For a finite V, all the remaining modes are bounded and the resulting mean field 
theory is well defined. The problem one is now faced with is the large system limit of the 
theory. It is easily seen from (34) that max^^o i| n |i^o a n = a(min ngZ D ,||n||^o \ \k n \\) diverges 
like V X I D as V — > +oo. This divergence jeopardizes the existence of the large system limit 
as it may cause an infrared divergence of V~ l 5^|| n ||<7/ (l ™! 2 )^- The problem is thus to 
determine whether the contribution of the infrared modes to V^ 1 ^2\\ n \\ <Vc {\a>n\ 2 )w remains 
finite when the discrete sum over n G 7h D is replaced by an integral over k G M. D according 
to (26). 

We consider the case in which Eq. (27) is replaced with 

,,(//)= { ! (35) 




where c > and k max > are fixed, and ((k) is a decreasing positive function of k, 
normalized to £(0) = 1, and such that ((k) > for k < k max and Cikmax) = 0. 
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For x > 0, Zi(x,k, (3V) in (29) is given by the first line of (31). In this case, the 
contribution of the infrared modes to (29) goes like / 5 k D ~ x dk if x > 0, and like k D ~ 3 dk if 
x = 0, where 0<(€l. On the other hand, for x < 0, the infrared contribution to (29), for 
which Zi(x, k, f3V) is now given by the second line of (31), behaves like \ \nk\~ 1 k D ~ 3 dk. 

According to these results, J-\(x < 0, /3V) diverges if D < 2. Thus, if a(k) behaves like 
1/k for small k, the large system limit of the mean field theory with e — 1 exists for D > 2 
only. 

Although our proof of existence of a stationary measure for the full dynamics (7) does 
not apply if a(k) is unbounded, it seems plausible that a stationary state exists for (34) 
if V is finite. Indeed, in this case the only unbounded mode is ao which is known to be 
modulationally unstable if its amplitude gets high enough. The larger |ao| the faster its 
modulational decay into bounded, smaller scale modes. Thus, it seems reasonable to expect 
that this transfer of energy toward smaller scales, and the subsequent dissipation according 
to the wave-breaking prescription, will eventually balance the growth of |a |, yielding a well 
defined stationary state. Proving this assertion comes within a study of the full dynamics 
(7) with wave-breaking prescription, which is the subject of a future work. 



V. SUMMARY AND PERSPECTIVES 



In this paper we have obtained stationary states for a system evolving according to a 
stochastic truncated NLSE on a D-dimensional torus, in the case of a focusing nonlinearity 
(A < 0). The stochastic dynamics are of a Langevin type (Ornstein-Uhlenbeck process) 
which act independently on every normal mode, both configurational and momentum, try- 
ing to bring its distribution to a Gibbs form at reciprocal temperature /3. Our approach 
consists in imposing a bound on the amplitude of each Fourier mode with either reflecting 
or absorbing boundary conditions (with reinjection at the zero amplitude point in the latter 
case). The latter boundary prescription is a simple modeling of a wave-breaking-like pro- 
cess, like the wave-breaking of Langmuir waves in plasmas. Since there is a finite number 
of modes, we are left with a system with a finite number of degrees of freedom evolving 
according to a stochastic diffusive dynamics on a compact phase space region. It follows 
readily that a stationary measure always exists, is unique, and is approached exponentially 
fast (Sec. III). We then studied the stationary state of a mean field version of (7). We have 
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shown that the mean field theory with A < admits a transition from a regime of low field 
values at small |A| (or large 0), to a regime of high field values at large |A| (or small (5) 
(Sec. IV). Field values at large |A| are significantly smaller with wave-breaking than with 
reflecting boundary conditions. 

The mean field theory fails to approximate the original model faithfully at the transition 
simply because, according to the result of Sec. Ill, coexistence of distinct stationary states 
is impossible in a finite volume for the stochastic model of Sec. II. This is of course similar 
to what happens with standard mean field theory in equilibrium where mean field predicts 
transitions where there is none in the real system [11]. Whether or not there is a counterpart 
of the transition for the full dynamics of Sec. II, as well as what the behavior of is like 
in the transition region (if it exists), are unresolved interesting questions that will probably 
require numerical simulations. 

Our modeling of wave-breaking is an oversimplified description of the actual process. 
For instance, in the wave-breaking of a Langmuir wave, not all the wave energy is transferred 
to the particles. This could be taken into account by replacing the delta functions on the 
right-hand side of (8) with smooth functions of small but finite width. We also make the 
particle degrees of freedom play the role of a reservoir for the waves. At first sight this seems 
reasonable since, due to the ultraviolet cut-off resulting from our wave-breaking prescription, 
the wave degrees of freedom are typically much fewer than the particle degrees of freedom 
[the ratio is of the order of (n e A^) _1 <C 1, where n e is the electron density]. Now, regarding 
the particles as a reservoir means that the particles should not be significantly affected by 
the violent exchanges of energy brought about by wave-breaking. In real situations, however, 
these exchanges generate a population of suprathermal electrons which can modify both the 
stochastic dynamics (7) and the wave-breaking itself. These intermittent processes must 
remain rare enough so that they have a chance to always occur in a (nearly) thermalized 
plasma. If this is not the case, then the particles cannot be idendified with a reservoir 
and (7) should be replaced with (i) a stochastic dynamics for the waves interacting with 
the particles, (ii) a stochastic dynamics for the particles interacting with the waves and a 
reservoir, and (iii) self-consistent wave-breaking prescriptions for (i) and (ii). 

Another important simplification is the fact that the wave-breaking for a given mode 
ignores the existence of all the other modes (it is just a copy of wave-breaking for a single 
monochromatic wave). This simplification is not very realistic as particles travel in wave 
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packets, rather than a monochromatic wave, in which several modes can have a high ampli- 
tude simultaneously. The wave-breaking of a given mode is therefore expected to be affected 
by the existence of the other modes via the perturbation of the electron motion they induce. 
What we need is a theory of wave-breaking for localized wave packets, or more generally for 
multimode structures. As far as we know this theory, which can be regarded as the highly 
nonlinear counterpart of linear Landau damping for localized wave packets, usually referred 
to as "transit-time damping" [12] [13], is still lacking [14]. Having such a theory would be 
of the utmost interest to improve our understanding of nonequilibrium stationary states for 
stochastic NLSE and other stochastic nonlinear wave equations of the same type. 
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